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ABSTRACT 

Oblique viewing sensors have been proposed for the Multlspectral Re- 
source Sampler (MRS), a follow— on satellite to the Thematic Mapper /Landsat D 
and proposed by the United States for launch in the late 1980' s. If an 
oblique viewing sensor Is to be Included on the MRS, then the potential in- 
formation in oblique measurements needs to be better understood. 

This paper analyzes the reflectance response with view angle of wheat, 
excluding atmospheric effects but otherwise simulating the response of a 
multlspectral scanner. The analysis is based upon spectra taken contin- 
uously in wavelength from 0.45 to 2.4 pm at more than 1200 view/ illumination 
directions using an Exotech model 20C spectral radiometer. Data were ac- 
quired six meters above four wheat canopies, each at a different growth 
stage . 

The analysis shows that the canopy reflective response is a pronounced 
function of Illumination angle, scanner view angle, and wavelength. The 
variation is greater at low solar elevations compared to high solar eleva- 
tions. 

INTRODUCTION 


Oblique viewing sensors have been proposed or scheduled for launch on 
at least two future earth resource satellite systems, the Systeme Probatolre 
d'Observation de la Terre' (SPOT) being developed in association between 
France, Sweden, and Belgium and the Multlspectral Resource Sampler^ (MRS) 
being developed by the United States. As these sensors are soon to be 
launched, the potential information in oblique measurements needs to be 
anticipated and better understood. 

The spectral flux sensed by these systems will be due to the absorption 
and bidirectional scattering of solar radiation by both the atmosphere and 
ground scene. Characterization and correction of the effect of the atmos- 
phere upon remotely sensed data have been considered elsewhere.' This paper 
analyzes the spectral bidirectional scattering properties of wheat with view 
direction, excluding atmospheric effects. 

If the ground scene is a level and a perfectly diffuse reflector, a 
Lambertian surface, then the response of the sensor will not vary with view 
angle, provided atmospheric effects may be neglected. However, never is an 
extended natural surface level and Lambertian; inevitable surface roughness 
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is responsible for light and dark surface areas and the consequent varia- 
tion in the sensor response with view angle. The roughness Inherent in a 
plant canopy provides an example of this process; the presence of shadows, 
cast by foliage, in differing proportions depending upon view angle sug- 
gests that the canopy is not a Lambertian reflector, that the sensor re- 
sponse will not be constant with canopy view angle. 

The reflective response of a vegetative canopy has been modeled*"** 
using physically based parameters such as leaf area index, leaf areas pro- 
jected In particular directions, probability density functions of leaf 
area with angle, probability of gap, reflectance and transmittance of can- 
opy components, index of refraction, etc. Analyses and parameter studies 
have investigated the predictions and properties of the models and, based 
on the models, offered insight into the canopy reflectance process . 

Very limited field verification of the •s- dels has been reported. 

Measurements of an information cxass such as wheat or corn made by 
mult Ispectral scanners mounted in aircraft often show large variation with 
scan angle. For example, during the Com Blight Watch Experiment, the 
size of the variation was sufficient to require preprocessing to remove 
the effect.*” The response of the aircraft sensors, which typically oper- 
ate over angular variations of *40® about nad^r, include variations due 
not only to bidirectional scattering by the ground scene but also due to 
the atmosphere. Consequently, data from the sensors cannot be interpreted 
as completely indicative of the ground scene. 

Measurements of the bidirectional scattering properties at large view 
angles for various plant canopies have been reported.*’** None of these 
studies involved measurements made continuously in wavelength from 0.4 to 
2.4 VHP and at a variety of view angles and crop growth stages. Nor have 
these studies extensively investigated the reflective properties of wheat, 
a crop of global economic importance grown worldwide. 

MATERIALS AND MErriODS 


Data were acquired on wheat (Trltlcum aestlvum L .) on four dates 
during 1975 and 1976 at Wllllston, North Dakota, USA (Lat. 48®8', Long. 
103®4<»') in support of the Large Area Crop Inventory Experiment (LACIE).** 
On each date agronomic measurements were made to characterize the condi- 
tion of the wheat canopy (Table 1) . Meteorologlc data (Table 1) were 
ac(|ulred at the North Dakota Agricultural Experiment Station at Wllllston, 
located near the test sites. 

All spectral data were acquired continuously In wavelength 0.46 to 
2.4 yn, using an Exotech model 20C spectroradlometer** (Exotech, Inc., 
Gaithersburg^ Maryland) positioned 6 m above the soil and bolted to a pan 

head mounted to a boom supported by a truck (Figure 1) . On each of the 

four dates, the truck and associated Instrument van and electrical gener- 
ator were located amidst the wheat and sufficiently distant from field 
edges to uniformly fill for all view directions the 15® field of view of 
the spectroradlometer. Spectral dats and a photograph of the Instrument 

field of view were taken in each of 33 view directions, eight azimuths 

(the eight points of the compass) at four zenith angles (15®, 30 , 45 and 
60®) plus nadir. On each day data were acquired approximately from three 
hours before solar noon to four hours after solar noon. The spectral data, 
acquired as radiances, were subsequently calibrated to bidirectional 
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reflectance factors (BRF) using spectral radiance measurements taken 
periodically of a 1.1 m x 1.1 m barium sulfate (BaSO^) painted field stan- 
dard.’® The reflectance of the field standard was measured with reference^^ 
to pressed BaSO^, a laboratory standard with known reflective properties. 

Spectral data and corresponding wavelength information were simul- 
taneously recorded on analog frequency modulated magnetic tape and on a 
multiple channel strip chart recorder. A crystal frequency controlled 
power source maintained tape transport speed specifications. On play back 
in the laboratory, the analog signals were sampled, digiti^od^ and then 
converted to bidirectional reflectance factors using algorithms imple- 
mented on IBM 370/1A8 and IBM 3031 computers. 

The field of view photographs, taken with each spectra, were used to 
assess spectral data quality. A spectra was discarded for analysis pur- 
poses if the associated photograph indicated the field of view of the in- 
strument might Include the shadow cast by the spectroradiomater, boom, or 
truck or in some way did not properly represent the scene. 

Analysis of the bidirectional response of the canopy was performed on 
data at 48 wavelengths selected at 0.02 pm intervals from 0.44 to 1.0 pm, 
at 0.04 prn Intervals from 1.0 to 2.0 pm, and at 0.08 tim intervals from 2.0 
to 2.24 pm. The wavelengths were chosen to concentrate analysis efforts 
in the visible spectral region. At a particular wavelength the spectral 
resolution of the data is better than 2.6% of that wavelength. 

At each of the 48 wavelengths, a stepwise forward regression program 
was used in a two part sequence to select from a global set of possible 
terms the twenty which best explained the variation in the BRF data with 
time and view angle (Figure 2). The global set included the terms of tht' 
spherical h..rmonic series’^ through , powers of time through L®, and all 
the interaction terms. The number of terms in the global set was first 
reduced by analyzing regression results at 10 wavelengths strategically 
placed across the spectrum and selecting from the global set those terms 
explaining most of the variation in the data. Because analysis of resid- 
uals indicated additional terms were needed to correctly model large var- 
iations in tbo BHr' in the* v»'*.ib1o spectral recion at view angles near 
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Figure 2. Coordinate system 
used for analysis. 
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nHdir, three terras, e.irli (invinp a Gaussian response with view .inple, were 
added to the select set of terms. Next a program for the regression of 
the data at each wavelength chose the best twenty terns from the select 
set and determined the coefficients of the regression equations. 

Analysis of the residuals of the twenty term regressions revealed no 
significant pattern as a function of the following variables; time, view 
zenith, view azimuth, sun zenith, sun azimuth, predicted BRF, and measured 
BRF. The coefficient of variation (R^ < of each regression varied syste- 
matically by date and wavelength, ranging between 0.86 and 0.96 In the 
visible spectral region, 0.93 and 0.96 In the near-infrared and 0.76 and 
0.96 In the middle- Inf rared . The standard deviation of the quantity (lOOZ 
residual/measured) varied between 3.5Z and 10. 7Z in the visible spectral 
region, 2.5Z and 5.2% In the near-infrared, and 3.3% and 16.4% In the 
middle- In fra red. 

The coefficients of the regression equations at the 48 wavelengths 
and four dates were processed by computer Into a computer program for 
exercizing the equations to perform parameter studies. Additional compu- 
ter programs were written to plot Che calculated values on a computer line 
printer. The resultant curves were traced and labeled by hand, photo- 
graphically reduced, and assembled Into arrays of plots. 

RESULTS 


The nonhalized BRF of wheat was plotted as a function of wavelength 
and view direction, at three times during the day, and for four crop 
development stages. The analysis assumes there are no atmospheric effects. 

Figure 3 shows the normalized BRF plotted as a function of wavelength 
for the wheat canopy measured 21 June 1976 three hours before noon looking 
90° to the sun azimuth. The BRF Is normalized to the BRF at nadir at each 
wavelength. The scale on the wavelength axis changes at 1.0 and 2.0 pm. 
Figure 3 shows that the BRF is a pronounced function of view zenith direc- 
tion and wavelength. In the green spectral region (0.56 pm) the BRF de- 
creases for view zenith directions between 0 and ab^ut 20° then Increaset 
for angles greater than 30°. In the red region (0.66 um) the BRF decreases 
until a view zenith direction of 40°, then Is constant. In the near- 
infrared (0.76-1.28 pm) the BRF increases with increasing view zenith 
angle. In the middle-infrared at 1.68 pm the BRF response Is like that in 
the green region while at 2.16 pm the response is similar to that In the 
red spectral region. 

The results shown in Figure 4 represent the normalized BRF of the 
canopy on one date 17 July 1976 three hours before noon. The results in 
Figure 4 are based on the simulated response of a conical scan, normalized 
to 1.0 at the canopy hot spot, the antlsolar point, and scanning in azi- 
muth angle with the zenith view angle fixed at the angle of the canopy 
hot spot. The results were obtained by dividing the wavelength axis Into 
three wavelength regions, visible 0.44-0.68 pm, near- Inf rared 0.76-1.28 pm, 
and middle-infrared 1.48-1.76 pm and 2.08-2.24 pm. On Figure 4 the range 
of the canopy BRF responses in each wavelength region is signified by the 
appropriate stippled area. For example, at an azimuth angle of -90°, the 
BRF In the visible region ranged between 0.5 and 0.55 of the BRF at the 
hot spot. Figure 4 shows that the response as a function of azimuth angle 
measured from the hot spot is symmetric and decreases. The decrease of 
the response is greatest In the visible region and least In the near- 
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infrared region. 


Figure S shows the wavelength and scan angle scales which applv to each 
plot In the 3x4 array of plots In Figures 6a and 6b. Tin w.ivelength scale 
changes at 1.0 and 20 piTu The scan angle scale Is organized for data from 
a nultlspectral line scanner scanning tbO® about nadir. Negative scan 
angles are zenith angle.s to the left of nadir. 

Figure 6 shows the BRF response of a MSS normalized to the response 
at nadir and scanning 160® about nadir. The wavelength and zenith scan 
angle scales of each of the 24 plots In Figure 6 are Illustrated In Figure 
S. Each line of plots In Figure 6 represents a particular time. Each 
column of plots in Figure 6 represents a specific crop development stage. 

In each plot the normalized BRF at a particular scan angle and wavelength 
Is Indicated by the height quantized in topographic notation by contour 
lines at 0.1 unit Intervals. The BRF on a contour line labeled 1.1 is 1.1 
times the BRF at nadir. The contour line Innnedlately adjacent to the line 
labeled 1.1 indicates a BRF value of either 1.0, 1.1, or 1.2 tiroes the 
nadir BRF, depending on the sequence of contour lines. 

The results in Figure 6a, the case of the MSS-equlpped aircraft fly- 
ing toward the solar azimuth, show the normalized BRF with scan angle is 
generally symmetrical about the nadir scan angle for all wavelengths, all 
tlirce times, and all four growth stages. For example, the BRF for [scan 
angle • t60°, ^3 hours, 21 June 1976, X-0.64 pm] Is 0.6 of the BRF at 
nadir at 0.64 uiri; similarly, the BRF at 0.8 piti Is 1.4 of the BRF at nadir 
for -3 hours and 1.5 of the BRF at nadir for +3 hours. In the near-infra- 
red spectral region, the BRF generally increases with increasing scan 
angle for all three times and all four crop development stages. In the 
visible and middle-infrared regions, a simple pattern doesn't exist: the 
BRF may increase or decrease depending upon wavelength, illumination angle, 
and development stage. 

Figure 6a also shows the BRF with scan angle for three hours before 
noon is very similar to the response three hours after noon, suggesting 
that the BRF is lalrly symmetric not only in scan angle but also 
with Illumination angles from noon. Generally, for a particular wave- 
length, whatever changes that do occur with scan angle are enhanced at the 
two times away from noon as compared to noon. In each plot the transition 
between the visible and near-infrared spectral regions Is abrupt for the 
green, healthy canopies (columns labeled 21 June 1976, 20 July 1975, and 
17 July 1976) and markedly less so for the senescent canopy (column la- 
beled 31 July 1976) without the strong red chlorophyll absorption band. 

Figure 6b, the case of the aircraft flying in a direction 90° to the 
sun azimuth direction, reve.ils a BRF similar to that of Figure 6a. The 
same transition phenomenon Is evident at 0.7 pin between the visible and 
the near Infrared regions. The curves are fairly symmetric with Illumina- 
tion angle from noon. In the near-infrared region, the BRF generally In- 
creases for increasing scan angles regardless of crop development stage 
or time. (Note the minor but systematic exception to this rule near 
nadir where the BRF Is between 0.9 and 1.0 of nadir BRF.) However, unlike 
Figure 6a the response curves of Figure 6b are not syiimietrlc with scan 
angle. 
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The results. Figure 4, arc consistent with the c>inccpt that the 
effects of shadowing on BRF are modulated by the effects of fight multiply 
scattered by canopy components. Radiation Is multiply scattered when It 
is reflected or transmitted more than cme t lm<> by foliage or soil fn the 
lanopy. Fight tends to be multiply scattered In crop canopies in regions 
ol the spectrum where the loliage absorbs little light (l.e., the near 
infrared band from 0.8 um to 1.1 iim) and tends to not be multiply scattered 
in regions of the spectrum where foliage absorbs a significant proportion 
of the Incident light (l.e., the red wavelength band prior to senescence). 
Multiply scattered light tends to reduce the contrast between two adjacent 
surfaces in the canopy, one surface Illuminated by direct solar radiation 
and the other surface shadowed. The BRF of a crop canopy Is the sum of 
the individual contributions of the shadowed and illuminated surface 
areas in the canopy. If the contrast between the shadowed and Illuminated 
areas is regllgible, then changes in the proportion of shadowed surface 
area to Illuminated surface area In the canopy will not be evident In the 
canopy BRF. Conversely, if the contrast is significant, then changes in 
the proportion of .shadowed to Illuminated surface area in a canopy will be 
evident a.s changes in the canopy BRF. Therefore, the canopy BRF is a 
function of both the proportion of shadowed to Illuminated surface area 
in the canopy and the importance to the canopy radiation environment of 
multiply scattered light, which reduces the contrast between the shadowed 
and illuminated surface areas. 

The results. Figure 4, are for one zenith view angle, that of the hot 
spot. If the canopy foliage were randomly distributed both spatlnlly and 
azimuthally, then, because the zenith view angle is constant, the relative 
proportions of the canopy components in the field of view of a sensor 
would not change with azimuth scar angle. The proportion of illuminated 
t(' shadowed foliage would be greatest at thr hot spot, an azimuth scan 
angle of zero, and least at an azimuth scan angle of 180 . If the BRF of 
the canopy were merely a function of the proportion of illuminated to shad- 
owed foliage, then the stippled areas of Figure 4 should coincide. That 
the regions do not coincide indicates other factors must be included In 
the analysis. 

At any azimuth scan angle away from the hot spot, the stippled areas 
are ordered from top to bottom, near-infrared, middle- inf rared , and visible. 
The ranking corresponds with a ranking of the importance of multiply scat- 
tered light in each spectral region. F.ach stippled area representing a 
spectral region is almost symmetric with azimuth scan angle from the hot 
spot; the proportion of illuminated to shadowed areas is similarly syttSBetrlc 
provided the canopy foliage is randomly distributed. Consequently, the 
results. Figure 4, are consistent with the argument that the effects of 
multiply scattered light serve to modulate the effects of shadowing on the 
BRF of the canopy. 

The results. Figure 3 and 6, ere consistent with the argument that 
shadows cast by foliage are an Important factor in the variation of BRF 
with view direction. Considering Figure 6a, the proportion of shadowed to 
sunlit foliage and soli should be approximately the same at equal angles 
left and right of nadir provided the foliage is randomly distributed. 

Sln.e the BRF of a crop canopy Is the sum of the individual contributions 
of the siiadowfd and illuminated surface areas in the canopy, then the BRF 
should be symmetrical in scan angle from nadir — as it is in Figure 6a, 
regardless of wavelength. Considering Figure 6b, the curves are asynmetrlc 
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In »c«n anglt about nadir. At all wavalan|thB, the normal lied aenaor ra- 
aponae la noticeably larger on the aide of the fllghtllne away from the sun, 
toward the canopy hot apot where ehndcwe caat by foliage would not be acen. 
Shadows cast by foliage would be observable on the opposite aide of the 
fllghtllne. the side with the lower reaponae Indicated by Figure 6b. Thua, 
the results support the argument that shadowing la an important factor In 
the variation of BRF with view direction. 

The results. Figure 6, are conalatent with the argument that ahadowa 
cast by foliage are an important factor In the variation of BUT not only 
with view direction but alao with Illumination angle. The proportion of 
ahadowed to Illuminated surface areas In a canopy changes with aun angle 
during the day, Tha proportion varies from unity (all ahadow) at eunrlac 
and sunaer to a low value aometlM during the day. For a particular view 
fixed relative to the Illumination direction, the proportionality, a func- 
tion in tlaw, will be ayamatrlc about aolar noon If the aclmuthal orienta- 
tion and spatial distribution of the foliage is eymmetrlc about a north- 
south line and If c"'’’? oundlng factors (l.e., photetroplea, plant geometry 
changes due to wlrd or moisture atreaa, etc.) are not Important. For a 
particular view direction fixed relative to the aun direction, the BIF 
should. In general, be aymnetrlc in time and sun angles about aolar noon, 
provided (1) the geometry of the canopy remains properly aymmetric through- 
out the day tnd (2) the apectrai properties of the canopy components are 
constant or vary symmetrically in time about aolar noon. (For example. If 
algnlficant soil surface dry down occurred during a day, then the canopy 
would fail tj satisfy criteria 2 when the soil reflectance changes. In any, 
were asynsetrlc about solar noon.) The results. Figures 3 and 6, show the 
BRF for a particuJ.ar wavelength and scan j ngle Is generally ayimsetric In 
time about solar noon, supporting the argument that shadowing la an impor-' 
tant factor for explaining the variation in BRF with Illumination angles. 

The results. Figures 3 and 6, are consistent with the Idea that the 
BRF changes in certain waya with aesn angle bacauae the probability of 
observing the varloua components of the canopy changes significantly with 
scan angle. For example, the probability of observing bare aoll la gen- 
erally greatest at nadir and decreases rapidly with increasing scan angle 
across the canopy. ' igure j and several plots In Figure 6 reveal an 
abrupt change of BRF w^.h scan angle near nadir for the visible spectral 
region, suggesting the effect of the aoll ahould be considered for under- 
standing the properties of the BRF under theae conditions. As a second 
example, for most canoplee at large scan angles only the upper layers of 
the canopy are visible. The upper layers of a canopy are well Illuminated 
— no higher foliage shadows the topmost layer. Thus, at large. Incremen- 
tally increasing scan angles the proportion of shadows should decrease and 
the BRF should Increase. Figures 3 and 6 reveal such a pattern In the 
near-infrared spectral region of all plots and In the visible portions of 
several plots. The pattern may exist throughout the visible spectral re- 
gion where the BRF might Increase at scan angles larger than the 60° 
measured. 
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Table 1. Ancillary meteorologic and agronomic data. 


Date 


Variable 

21 Jun 76 

20 Jul 75 

17 Jul 76 

31 Jul 76 

relative humidity (%) 

51 

54 

36 

57 

air temperature (C°) 

19 

27 

28 

23 

barometric pressure (can Hg) 

770.9 

759.4 

771.4 

777.2 

cloud cover (2) 

1 


1 

5 

wind direction 

northeast 

southwest 

southeast 

southeast 

wind speed (km/hr) 

14 

16 

10 

13 

cultlvar 

Waldron 

Wells 

Ellar 

Ellar 

maturity stage* 

3.5/boot 

4.5/fully 

5.1/milk 

5.4/ripe 

row dlrectlor 

east-west 

headed 

east-west 

north-south 

north- south 

row width (m) 

0.18 

0.21 

0.18 

0.18 

fruit count (p'jr m^) 

0.0 


444.4 

394.4 

plant count (per m^) 

477.8 

310.0 

455.6 

405.6 

plant heigh (m) 

0.48 

0.72 

0.85 

0.86 

leaves per plant 

5.0 


4.0 

4.0 

leaf condition (2) 
green 

93 


27 

0 

yellow 

3 


7 

0 

brown 

4 


66 

100 

dry biomass-total (gr/m^) 

216.9 

345.1 

689.8 

625.9 

fruit 

0.0 


268.1 

330.2 

green leaves 

84.7 


51 

0.0 

yellow leaves 
brown leaves 
stems 

132.3 


25.4 

345.3 

56.8 

238.9 

fresh biomass-total (gr/m^) 

1131.1 


1466.1 

840.0 

plant molstute (X) 

81 


53 

25 

leaf area index** 

1.85 

1.48 

0.81 

0.0 

♦maturity stage according to 

Large** 




**leaf area index is the green 

one-sided leaf area per unit ground area 



